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INTRODUCTION 


During the final free flight (FF-5) of the shuttle's approach and landing test 
(ALT) phase , the vehicle underwent pilot-induced oscillations (PIO's) near touch- 
down (refs , 1 to 3) . The oscillations were present in both the pitch and roll axes 
and were initiated when the pilot made pitch controller inputs in an effort to control 
sink rate by changing pitch attitude. Because the control inputs were large and 
fairly rapid , the elevens rate limited in the pitch axis at the maximum priority rate 
limit set in the computers. The eleven rate limit also limits the vehicle's roll control 
capability , and this was partially responsible for the lateral control problem . 

Several unpublished studies indicate that time delays as well as priority rate 
limiting were a significant factor in the PIO's. A simulator study of the effect of 
time delays on shuttle PIO's is reported in reference 4. 

This report describes the combined effect of pilot input rate limiting and time 
delays . Frequency responses are predicted for various parameters under rate 
saturated conditions by using nonlinear analysis . 


SYMBOLS AND ABBREVIATIONS 


Physical quantities in this report are given in the International System of Units 
(SI) and parenthetically in U .S . Customary Units . The measurements were taken 
in U.S . Customary Units. Factors relating the two systems are presented in 
reference 5 . 
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moment of inertia about y-axis , N-m (slug-ft ) 

imaginary part of root , per sec 
gain constant 

gain factor out of the shaping function 

Mach number 
pilot-induced oscillation 
pitch rate , rad/ sec 

2 2 

dynamic pressure , N/m (Ib/ft ) 

reference area, m^ (ft^) 

Laplace transform variable , per sec 
time delay , sec 
transfer function 
time , sec 

velocity , m/ sec (ft/ sec) 
weight, kg (lb) 
angle of attack , deg 
transfer function denominator 
elevon deflection , deg or rad 

commanded input , rad 
pilot input , deg or rad 
shaping function output , deg or rad 
trim attitude , deg or rad 
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0 

pitch attitude , deg or rad 

T 

time constant, sec 

9 

phase angle , deg 

CO 

frequency, rad/ sec 

Subscripts: 


CP 

cockpit 

c 

commanded 

e 

error signal 

P 

pilot 

R 

generalized input 

ref 

reference 

1,2,3,... 

sequential calculation 


A dot over a quantity denotes the first derivative with respect to time . 

VEHICLE DESCRIPTION 


Aircraft and Test Configuration 

The space shuttle vehicles were built by Rockwell International under a 
contract awarded by the National Aeronautics and Space Administration. The first 
vehicle, which was designed for unpowered flight, was to be flown in approach and 
landing tests only . 

The shuttle used for the approach and landing tests was carried aloft on the 
fuselage of a modified B-747 airplane and then launched at an altitude and position 
within gliding distance of a dry lakebed at Edwards Air Force Base. Five glide 
flights were made with this vehicle . A special tail cone was installed to cover and 
streamline the aft portion of the vehicle . The first three flights were flown with 
the tail cone installed . On the last two flights , the tail cone was removed and 
dummy engines were installed to simulate the drag or aft profile that would exist 
on the rocket-launched shuttle during the orbital flight tests . 

The ALT configuration has essentially the same flight characteristics as the 
orbital flight test configuration . There are slight differences in the aircraft 
dynamics due to the differences in weight and inertia and the thermal protection 
system on the reentry vehicle . This study is limited to the ALT vehicle . 
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Figure 1 is a three-view drawing illustrating the relative size and location of 
the primary control surfaces of the shuttle . Pertinent physical characteristics are 
presented in table 1 . 


Primary Flight Control Systems 

Aerodynamic surfaces . — Full span independently actuated elevens provide both 
pitch and roll control. A conventional rudder mounted on the trailing edge of a 
single vertical stabilizer provides directional control. Speed brakes are activated 
by mechanically splitting the rudder symmetrically . A flap on the lower aft fuselage 
provides active trim control. 

Control systems .— A digital fly-by-wire control system is implemented on the 
space shuttle. All elements of the control system , including the hand controller 
electrical signals , sensors , computers , and electrical paths up to the servoactuators , 
are quadruply redundant. Three independent hydraulic systems (APU’s) supply 
power to the surface actuators. The independent operation of the four control systems 
and three hydraulic systems provides fail-operational/fail-safe capability . 

Rotational hand controller . — Both the commander's and the pilot's station are 
provided with identical rotational hand controllers. In pitch, the controller pivots 
about the palm of the hand; in roll, it pivots about a point slightly below the base. 
Therefore , the pilot effectively applies a torque gradient rather than a force gradient . 
The commander's and the pilot's station controllers are not mechanically linked, but 
if both are deflected the input signals are combined and the elevens are deflected 
accordingly . However , the controllers are spring centered and null at zero 
electrical output . 

Pilot control and augmentation .— A pitch rate command system is used on the 
space shuttle for the approach and landing phases of flight . The system is imple- 
mented as a digital flight control system resident in the onboard flight control 
computers . The mathematical block diagram in figure 2 shows the major elements 
of the closed-loop pitch attitude control system . The rotational hand controller input , 
6 , commands an output given by the pitch stick shaping function . The shaping 

function output ,6 , is prefiltered to form the pitch rate command , 0 . In the for- 

0 c 

s 

ward path and downstream of the error signal, a lead-lag filter along with a gain 
term, f(q) , gives the desired pitch rate response. The gain is scheduled as a 
function of dynamic pressure (200/q =1.64 for touchdown) . Pitch rate command is 
achieved by positive feedback around the power actuator and servo, yielding an 
equivalent system for the actuator loop with a "free" s in the denominator . As such , 
the error signal , which is made up of the commanded response and aircraft pitch 
rate , is continually driven to zero . Therefore , the integration rate is proportional 
to the scheduled gain. Internal to the actuator and servo loop, the priority rate at 
which the actuator can move is restricted to 20 degrees per second , and the travel 
limit, which has been changed, is restricted for the following calculations to 21.5®. 
The filters , gains , summations , restrictions , and pitch stick shaping function are 
implemented in the shuttle's digital flight control computers. 
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A more complete description of the total flight control system and its design is 
given in reference 6 . 


ANALYSIS AND DISCUSSION 


Shuttle PIO Experience 

The first four landings of the ALT shuttle configuration were made on the lake- 
bed without incident , as reported in references 2 and 3 . The handling qualities 
were judged to be good; longitudinal control was rated approximately 2 on the 
Cooper-Harper scale (ref. 3) . 

A precision landing on the runway was attempted during the last flight. During 
the landing , the pilot induced both a longitudinal and a lateral PIO . The PIO first 
occurred in the pitch axis and then propagated to the roll axis . The oscillation was 
initiated by pitch controller inputs that were made to control sink rate . The exces- 
sive controller inputs caused the elevens to rate saturate . A right roll disturbance 
was then introduced as the main gear touched the runway . A lateral PIO then 
followed because roll control was restricted by priority rate limiting and control 
authority , which favored the pitch axis . The combined time of the oscillations was 
more than 12 seconds , and both oscillations were near the same frequency . Because 
the oscillations originated in pitch, it was assumed that the PIO problem was pri- 
marily a longitudinal control problem , with the roll control problem system induced . 
The flight conditions near the time of the PIO's are listed in table 2. 

Figure 3 presents a time history illustrating the longitudinal PIO . Pilot inputs 
to the hand controller, 6 , were approximately 0.5 hertz and in excess of 20° peak 

to peak. The average travel of the elevens was about 15° peak to peak. In addition, 
the control surface time history indicates that the actuators had reached the priority 
rate limits . The resulting pitch rate during the longitudinal PIO reached a peak to 
peak value of 7 degrees per second and appeared to be neutrally stable . 

Prior investigations had predicted the shuttle to be susceptible to PIO's. An 
analytical study based on the shuttle's aerodynamics concluded that the shuttle 
would be prone to a longitudinal PIO (ref. 1) . Also, simulator studies conducted 
concurrently with the approach and landing tests , as reported in reference 2 , 
showed a tendency towards PIO's. 

After the PIO's were experienced during FF-5 , the rotational hand controller 
signal shaping was modified to include a suppression filter like that reported in 
references 7 and 8 , and an analytical study was undertaken to investigate and 
evaluate the following: 

1 . System lag caused by rate limiting in the actuator loop . 

2 . The effect of varying the rotational hand controller input on the total atten- 
uation and phase lag . 

3 . The effect of time delay in the augmentation loop on the gain and phase lag . 
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The following section gives an analytical description of the long^itudinal or pitch 
control mode implemented on the shuttle . 


Pitch Control Modeling 

Presented in figure 4 is a conceptual block diagram of the shuttle control system . 
The pilot attempts to minimize pitch attitude error; he also tries to maintain the 
desired sink rate by changing the vehicle’s reference attitude. The amount of 
rotation or torque applied by the pilot to the controller per unit error signal is 
referred to as pilot gain. 

Figure 5 shows the parabolic hand controller shaping function programmed on 
the shuttle . The primary purpose of the shaping function , as it is in most fly-by- 
wire systems , is to increase the gain as a nonlinear function of controller displacement . 
Figure 6 presents the gain variation as a function of controller input , 6^ . In the 

P 

command path and downstream of the shaping function, a gain constant of 0.4 
degrees per second per degree is scheduled for the final approach and landing 

A6 

e 

g 

phases of flight. The pitch rate commanded, therefore, is 6^ — 0.4. To 

P 

^ P 

improve the handling qualities , a lead-lag filter ,1,5 /g + 2*7) ’ cascaded in the 
forward path of the augmentation loop . 

It was reasoned that dead time resulting from transport lag plus the delay due 
to the rather low rotational hand controller sample rate of 12.5 cycles per second 
had worsened the PIO problem. Consequently, in the analysis a time delay function 

e was added into the forward path of the augmentation loop . To assess the effect 
of this type of nonlinearity on the closed- and open-loop behavior of the system , the 
time delay function, T, was assigned various values: 0, 0.1, 0.2, and 0.4 second. 
Pitch rate command was achieved by positive feedback around the power actuator 
and servomechanism , as shown in figure 7 . The system for the actuator loop with 
no rate limiting is given by 


_ (s/1.5 + 1) 

8 s(s/30 + 1) 

e 


Therefore , the integration rate was proportional to the 


f(q) = 


200 

q 


1.64. 


gain f (q) (fig . 4) , where 


During FF-5 , the longitudinal control system underwent an extensive amount of 
rate limiting prior to touchdown , as is evident in the elevator time history of figure 3 . 
The boundary limits of both position and rate limit as a function of input amplitude and 
frequency are presented in figure 8. Three 5 inputs are shown parametrically as 

% 
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a function of frequency . These variations were obtained from a linear calculation 
with the shaping gain equivalent to the peak controller input . At realistic trim 
conditions and from controller inputs, 8 , up to 15°, the system would only 

% 

experience priority rate limiting over the frequency range shown. For controller 
inputs typical of the PIO condition of FF-5 (6 = 10°) , the system would rate limit 

% 

at 2 . 5 to 3 . 0 radians per second . Because rate limiting would be expected to add 
phase lag into the overall system , the actuator loop was first investigated in a 
separate or isolated fashion. 

Presented in figure 9 is a time history typical of the digital simulation of just 
the actuator loop . For this particular example , the input amplitude to the actuator 
was 10° at a frequency of 3 radians per second. The system rate limited at 20 degrees 
per second over 90 percent of the time , and the through gain based on peak values 
was less than unity. Several computations were made at various frequencies and 
amplitudes to obtain a broad comparison to a model of the actuator loop . 

Figure 10 compares both the amplitude ratio and the phase angle of the digital 
simulation to the algorithm that describes the actuator model . For a value of 6 

^c 

of 10° (near the PIO condition) , the rate limit causes attenuation from 3 radians per 
second on . Additional phase lag is also added into the system from there on . The 
' algorithm describing the actuator model is developed and presented as a subroutine 
in appendix A . 


Integrated Analysis of the PIO 

Nonlinear method . — During the PIO , an extensive amount of rate saturation 
prevailed. In addition, it was estimated (as reported in ref, 2) that up to 0,38 sec- 
ond of equivalent time delay may have been present in the augmentation loop (fig. 4) 
Therefore , a realistic analysis requires consideration of these two nonlinearities as 
well as the gain variation through the shaping function . The following discussion 
utilizes a nonlinear computation of the system gain and phase as a function of pilot 
input and frequency . 

The method , in general , is simply to assign a value to the error vector , which 
has amplitude and frequency but phase equal to zero. Then an input/output vector 
can be computed successively through each element , keeping track of amplitude and 
phase algebraically until an open-loop vector for the augmentation loop is obtained . 
The input to the augmentation loop is then determined by adding the feedback vector 
to the error vector . The reciprocal or inverse functions in the command path 
(reversing the direction of flow) are used to obtain the pilot input, 6 . This then 

P 

gives a closed-loop vector relationship of the desired elements as a function of 6^ 
and frequency . p 

Appendix B outlines the nonlinear computational procedure . The aircraft 
aerodynamics and transfer functions are presented in table 3 . 
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Elevon attenuation and system response characteristics .— If the amplitude of 
the input is great enough, elevon travel will be attenuated due to rate limiting as 
frequency is increased. Figure 11 presents the elevon travel as a function of con- 
troller input , 6 , and frequency . Also shown are the boundary curves for minimum 

G 

p 

perceptible and fully developed rate limiting. For inputs of 5°, no attenuation is 
expected because of rate limiting for frequencies up to 8 radians per second , as 
indicated by the circle symbols , which remain below the minimum rate limit boundary . 
At 8 = 10° the elevens are fully limited in a sawtooth fashion at frequencies from 3.5 

% 

radians per second on. If the input increases to 15° , the rate limiting begins at 
approximately 2 radians per second and then follows the fully developed rate limit 
boundary identical to the 10° input. 

The effect of time delay on elevon amplitude as a function of frequency is shown 
in figure 12. Time delay is varied parametrically from 0 to 0.4 second. A constant 
controller input amplitude of 10° is applied to the shaping function. Increasing the 
time delay increases the elevon amplitude as frequency increases . From 4 radians 
on , rate limit tends to limit the elevon amplitude , regardless of the amount of dead time . 

The normal force response of the aircraft , particularly the crossover phase 
frequency , is believed to affect pilot ratings . An assessment of the PIO tendencies 
based on both pitch attitude and the aircraft's normal force response is reported in 
reference 1 . Presented in figure 13 is the effect of time delay on the frequency 
response of a^ at the cockpit per controller input, A value of 5° was used for the 

controller input in order to stay below any rate limiting over the frequency range 
shown . It is clearly evident by the increase in gain and phase lag that an excessive 
amount of dead time (T = 0.4 sec) causes the aircraft to be more responsive per 
controller input. The peak increase occurs at frequencies from 2.5 to 3.0 radians 
per second. For a time delay of 0.4 second the crossover phase frequency (fig. 13(b)) 
is about 2 radians per second; with no time delay , the crossover phase frequency is 
above 10 radians per second. The same parameters are shown in figure 14 for three 
controller inputs at a constant value of dead time , T , of 0 . 1 second. For a controller 
input of 5° , the system behaves in a linear fashion. At 10°, the response is attenu- 
ated by rate limiting from approximately 2 . 5 radians per second on . As shown in 
figure 14(b) , rate limiting also adds phase lag into the system. 

The calculated frequency response of the shuttle's pitch rate command system 
is presented in figure 15. A constant value of 8 of 5° is applied to the shaping 

function. The effect of dead time is shown parametrically and is varied from 0 to 0.4 
second. Again, and consistent with normal force response characteristics, the 
system would be highly susceptible to oscillations at approximately 3 radians per 
second with 0.4 second of dead time in the augmentation. Figure 16 shows the 
frequency response of the pitch rate command system for a constant value of dead 
time of 0 . 1 second at three values of pilot input . For a value of 8 of 5° , the 

P 

system retains linear behavior throughout the frequency range shown. For values 
of 8 of 10° and 15°, rate saturation limits the attenuation, as shown by the diagonal 

% 

lines. The breakpoint for a value of 8 of 10° is approximately 3.5 radians per 

®P 
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second; it is approximately 2 radians per second for 8 = 15°. Rate saturation 

% 

adds about 50° of lag into the system, as shown in figure 16(b) . 

Figure 17 presents the computed results of pitch attitude response per pilot 
input versus frequency. Three values of controller input are presented. As pointed 
out previously , the system behaves in a linear fashion for 8 = 5° . The amplitude 

% 

rates normally decrease at 20 decibels per decade until the system starts to rate limit. 
The breakpoint occurs at 2.0 and 3.5 radians per second for values of 8 of 15° 

% 

and 10°, respectively. The crossover phase frequency is approximately 3.5 radians 
per second for both 5° and 10° of 8^ (fig. 17(b)) . About 70° of additional phase 

P 

lag is added into the attitude system at maximum pilot inputs . 

In order to sustain an oscillation in a closed-loop system , the total gain must be 
unity at a phase shift of -180° (Nyquist criteria) . Assuming that the pilot acts as 
pure gain and with no phase changes, the pilot gain required to cause a PIO would 
be the gain margin at a -180° phase shift. Because of system nonlinearities, the 
critical pilot gain is also a function of controller input , 8 . Figure 18 presents 

PIO frequency and critical pilot gain variation as a function of pilot input for various 
values of dead time . Increasing the pilot input decreases the pilot gain required to 
PIO . This effect is due primarily to the increase in gain across the shaping function 
with an increase in 8^ (fig. 6) . Increasing the amount of dead time in the augmen- 

P 

tation loop reduces both the crossover frequency and pilot gain required to PIO . For 
reasonable values of dead time (T = 0 . 1 sec) , the critical pilot gain would be reduced 
at least 50 percent. 


CONCLUSIONS 


A longitudinal PIO experienced on the shuttle was analyzed by a nonlinear tech- 
nique . In general , the analysis indicates that: 

1 . For controller inputs of 5° or less with no dead time , the control system 
behaves in a linear fashion. 

2 . For 10° of controller inputs , regardless of the amount of dead time , the 
eleven loop rate saturates at frequencies above 4 radians per second. 

3. The pitch attitude crossover phase frequency is about 3.5 radians per second 
for both the 5° and 10° controller input. 

4. At the PIO condition, rate limiting in the forward path decreases the system 
gain and also adds phase lag into the system . 

5 . Increasing the amount of dead time reduces the phase crossover frequency 
and reduces the pilot gain required to PIO . An excessive amount of dead time 
(0.4 second, for example) significantly increases the PIO tendencies. 
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APPENDIX A- MATHEMATICAL DESCRIPTION OF THE ACTUATOR LOOP 


A pitch rate command system was used on the shuttle for the approach and 
landing phases of flight . The system was designed such that the error signal in 
the augmentation loop was continuously driven to zero so that the output or vehicle 
pitch rate actively tried to match the commanded input (type 1 system) . To effect 
the necessary integration , positive feedback was implemented around just the actu- 
ator loop , as shown in figure 19 . To provide lead in the system , a first order lag 
was programmed in the actuator loop feedback path. Also, the preselected or 
priority rate limit plus the travel limitations were programmed in the forward path 
of the actuator loop . However , for this particular investigation , travel restriction 
was not necessary; consequently , for this study that gain factor was always unity . 
Because of the nonlinear nature that existed in just this part of the system , it was 
felt worthwhile and convenient to define and investigate the actuator loop separately . 

In general, what was desired was an input/output relationship that was a function 
of the same input . In the development of the nonlinear algorithm , it was found 
expedient and useful to begin with a linear input/output relationship . This was 
accomplished by iterating around the loop with both nonlinear elements , G6 and G7 , 
equal to unity . With reference to the following subroutine , consider first the input 
signal XX (line 7) , which is a vector having amplitude , phase , and frequency . 


SUBROUTINE TING{XX#Z) 

COMHOM /TFS/S 

COHNON /JAKE/X(5)f ZX>YNEU«Yf W#A#PHI»XX2 
COMPLEX S#XX#Y«Z*U»X#YNEW#G6»67 
3 COMPLEX ZX»XX2 

N*100 
Y-XX 

DO 10 K«1#M 
X(1)-Y+XX 

10 X(Z)« S^XCl) 

X(3)«X(2) 

X(^J- (1/S)*X(3) 

Xt5)« X(4) 

Z- C20/(S«20) )tX(5) 

15 YNEW*U«5000/(S<i>1.50000) )PZ 

1P(ABS(CABS(Y)-CA6S( YNEV) ) .LE.0.00001) GQ TO 11 

Y - YNEW 

10 CONTINUE 

11 CONTINUE 

20 XC2)- S«XC1) 

X(3)-G6(XC2)) 

X(^)- (1/S)»X(3) 

A«CABS(X(A) ) 

SIG -REAL(XC4) ) 

25 OMG -AIMAG(X(4) ) 

PHI - 1.570796 
IF(SZG.EQ.O.O) GO TO 6 
PHI •ATAN2(0MG#S1G) 

6 CONTINUE 

30 DEL - 0.0 

IF(CABSCX(1)1*AIMA6(S).LE..5^8) GO TO 7 
VV - .54a/<CABS(X{l) )«AlMAG(Sn 
DEL - 1.570796*VV -1.570796 

7 CONTINUE 

35 TPD «PHI ♦ DEL 

U - A* SIN(TPD» 

Y • A*C0SCTPD) 

X(4I* CMPLX (V#U) 

X(5)* XU) 

^0 2- (20/CS420) )«X(5) 

Y -Cl. 5000/(S*l. 50000) )*Z 
W • Z/CXCD-Y) 

Z -W*XX 
A-CABS(W) 

45 SXG-REALCU) 

OMG-AIMAG(W) 

PHI-90 

XF(SIG.EQ.O.O) GO TO 12 
PHI- ATAN2(0M6#SI6)*57.296 
50 12 CONTINUE 

ZX-Z 
XX2-XX 
RETURN 
END 
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Next , sum the input signal with Y and initialize for the first iteration by setting 
Y equal to XX. The vector X(l) will then be differentiated to obtain the rate vector 
necessary for subsequent computations . Throughout the first series of iterations 
the rate restriction G6 is equal to unity gain, or X(3) = X(2) . The vector X(3) is 
integrated to get back to a position vector X(4) . Since there is no appreciable 
amount of travel limitation, the vector X(5) will equal the vector X (4) , or the gain 
of G7 will be unity. A first-order lag was used to approximate the power actuator 
characteristics . That is , 


and 


A test follows at this point to determine whether the desired convergence 

_5 

(1 X 10 ) has been obtained. If not, Y is replaced by YNEW and the procedure is 

repeated until the desired convergence is reached . When the "if statement" is satis- 
fied , the computation continues or leads into the nonlinear computation (FORTRAN 
statement 11 on ) . The final value of the vector X(l) from the previous iteration is 
differentiated to get X(2) . The rate restriction is then applied to X(2) . In other 
words, X(3) = G6 (X(2)). 


However , the maximum gain for the sawtooth type of wave form rate limiting 
is used for these computations, as follows: 


X(3) < (20/57.3)71/2 = 0.548 rad/sec 


The rate limit (0 . 548) is entered into the following function subprogram . 


COMPLEX FUMCTIDN G6(X) 
COMPLEX X,CLIM 
COMPLEX A 
A«X 

A«CLIM(A,,5^6) 

66>A 

RETURN 

END 


A general limiting algorithm is represented by another function subprogram . 

COMPLEX FUNCTION CLIM(A#X) 

COMPLEX A>B 
Y-CABS(A) 

IFtY.EQ.0.0) GO TO 1 
B«A 

IF(Y.GT.X) 6«(X/Y)*A 
CLIM-B 
RETURN 
1 CONTINUE 

CLIM«CMPLX(X>0.0) 

RETURN 

END 
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Integrating X(3) gives a new position vector X(4) . Next, the absolute value , 
real, and imaginary parts plus the phase angle of X(4) are determined. From 
FORTRAN statement 6 on down to statement 7 , a new phase angle is computed based 
on the intersection of a sinusoidal oscillation and a rate limited wave train oscillation . 
This new angle (TPD) is then used to compute a new vector for the forward path. A 
final computation gives the vector amplitude ratio , W , representing the nonlinear 
system as a function of the initial input vector , XX . 
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APPENDIX B- BASIC NONLINEAR ANALYSIS AND COMPUTATIONAL PROCEDURE 


The flow diagram shown in figure 20 represents the shuttle's basic control 
system , aerodynamics , and normal acceleration at the cockpit . 

In the diagram, T1 represents the pilot transfer function, and for this study T1 
will be considered a pure gain term . T2 is the nonlinear shaping function shown in 
figure 5. All transfer functions or the inverse functions of Tl, T2, T3, and T4 are 
presented in tables 3 and 4, except for the actuator algorithm T7, and that function is 
derived and presented as a subroutine in appendix A . 

To obtain an open-loop frequency response for the pilot loop , an absolute veilue 
is assigned to XI , which will be the pilot input to the shaping function for unity 
pilot gain. The task in general is to find an absolute value for X2 that is close to 
XI. For example; let 6 equal 5° of pilot input, or XI equal 0.0873 radian. The 

desired accuracy is ABS(CABS(X1) - CABS(X2) < 1 X 10 

The first iteration is initiated through a "do loop" by assigning a number to the 
real part of X7 , such as 

1 = 1 , 1000 

SIG Z7 = 1.0 X I 

Next, the vector signal is traced around through each block, including the 
inverse functions to the left of the augmentation summation junction. 

Now , if and when CABS X2 > XI, the previous computation for I is used. So 
RBASE = I - 1 or a new index is J = 1, 10 and SIG X7 = RBASE + J(O.l) . 

The above scheme is repeated six times , moving to the left one decimal point 
at a time and proceeding with a new "do loop" until the desired accuracy is obtained. 
The complex value of each vector is successively computed through each block and 
summation junction. A final computation then yields the input/output amplitude and 
phase relationship as a function of pilot input and frequency . 

The logic , vector arithmetic , and transfer functions shown in the tables and 
appendix A for various configurations were programmed on a digital computer . In 
the equations , s = j to , where cd is a discrete frequency . 


Dryden Flight Research Facility 
Ames Research Center 

National Aeronautics and Space Administration 
Edwards, Calif. , July 29, 1981 
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TABLE 1. -PHYSICAL CHARACTERISTICS 


s, 

c, m (ft) 

«ref’ 

w , kg (Ib) 

eg, percent of . . 

Iy> kg-m^ (slug-ft^) . 


249.91 (2690) 

12.06 (39.56) 

32.77 (107.5) 

. . 68,038.50 (150,000) 
66.25 

6,833,703.42 (5,040,282) 


TABLE 2. -FLIGHT CONDITIONS NEAR TIME 
OF PILOT -INDUCED OSCILLATIONS 


a, deg 7 

M 0.32 

h , m (ft) . <694.94 (<2280) 

P 

q, N/m^ (Ib/ft^) . . 6654 to 5841 (139 to 122) 

V, m/sec (ft/sec) 97.54 (320) 

0Q, deg 7 


TABLE 3. -AERODYNAMIC STABILITY COEFFICIENTS 
AND TRANSFER FUNCTIONS 

Body axis stability coefficients: 


= 0.0053 
a 

C^ =0.0016 
6 

e 

C^ =0.048 

C, =0.018 
^6 

e 


C =-2.00 
m 

q 

C =0.001 
m 

a 


m. 


-0.0088 


Transfer functions: 


u _ 2.44 (s + 18.48) (s + 0,436 ± j 0.649)^ 


a _ -0.226 (s + 5.99) (s + 0,0108 ± j 0.0770)' 


q _ -1.30 s(s + 0.648) (s + 0.0349) 

6 A 

e 


n _ 2.26 (s + 2.089) (s - 1.754) (s + 0.0277) (s - 0.0156) 


where 

A = (s + 0.0945 ± j 0.122)^ (s + 0.887) (s - 0.100) 


TABLE 4. -TRANSFER FUNCTION DEFINITIONS 
(IXll designated input; (|X1| - |X2|) < K; 
see table 3 for aerodynamic transfer functions] 


Function 

Vector output 


XI = X2 (open loop) 

T1 (inverse; pilot 
transfer function) 

X2 = * X3 , T = 0, T = 0 

T2 (inverse hand 

A = 1X41 - , 

controller) 

B = (13.83 + 20.66 * 57.3 * A)”’'* 

C = -2.569 + B 

X3 = (C/A * 57.2) * X4 

T3 (inverse) 

yA _ S + 50 3k V- 
0.05 

T4 (inverse) 

X5 = §4-r4 * 0-667 * X6 

T5 

X8 = 1.64 * X7 

T6 

X9 = e'^*® * X8, T = 0.1, 0.2, 0.4 

T7 

Subroutine TING (XX, Z) (app. A) 

T8 

X11=1.0*X10 „ 


If ixill > 0.375, XU - * XU 

T9 

X12 = * XU 


o 

e 


a 

TIO 

X13 = ^*X11 


o 

e 

Til 

X14 = ^ * XU 


°e 

T12 

X15 = • XU 


^e 

T13 

X18 = s!°5o*X15 

T14 


T15 

X20 = - * X15 

S 

T16 

X16 = (1.55 * s) *X15 

T17 

X17 = 1.0 ♦ X13 


X6 = X19 - X7 \ 
X21 = X16 X171 




00 


23.79 

( 78 . 06 ) 





Figure 1. Three-view drawing of shuttle 
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Dimensions in meters (feet) 




SYSTEM LIMITS 
RATE = 20 ® /sec 
TRAVEL - 21 . 5 ® 



Figure 2. Mathematical description of shuttle pitch attitude control system. 
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Figure 4. Conceptual block diagram 
closing the pitch attitude loop. 
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of shuttle control system with pilot 




























£o, rad/sec 

Figure 11. Eleven amplitude attenuation due to rate limiting as a 
function of frequency and pilot input. T = 0.1 second. 
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CO , rad/sec 

Figure 12. Eleven amplitude attenuation due to rate limiting as 
a function of frequency and time delay. 6 = 10°. 
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CO, rad/sec 
(aj Amplitude ratio. 

Figure 13. Effect of time delay on normal acceleration at 

the cockpit. 6 = 5°. 
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per sec 



CO, rad/sec 
(a) Amplitude ratio. 


Figure 15. Effect of time delay on pitch rate 
command system. 6 =5°. 
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CO, rad/sec 
(a) Amplitude ratio. 


Figure 16. Effect of pilot input on pitch rate 
command system. T = 0.1 second. 









Figure 19. Actuator loop mathematical model 
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